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ABSTRACT 

We present spectropolarimetric observations of the FK Com star HD 199178 obtained 
between 1998 December and 2003 August at the Telescope Bernard Lyot (Observatoire 
du Pic du Midi, France). We report the detection of a photospheric magnetic field and 
reconstruct its distribution by means of Zeeman-Doppler Imaging. We observe large 
regions where the magnetic field is mainly azimuthal, suggesting that the dynamo 
processes generating the magnetic activity of HD 199178 may be active very close 
to the stellar surface. We investigate the rapid evolution of surface brightness and 
magnetic structures from a continuous monitoring of the star over several weeks in 2002 
and 2003. We report that significant changes occur in the distribution of cool spots and 
magnetic regions on typical timescales of the order of 2 weeks. Our spectropolarimetric 
observations also suggest that the surface of HD 199178 is sheared by differential 
rotation, with a difference in rotation rate between equatorial and polar regions of the 
order of 1.5 times that of the Sun. 

Key words: Line : polarization - Stars : rotation - imaging - activity - magnetic 
fields - Stars : individual : HD 199178. 



1 INTRODUCTION 

Our understanding of the Solar dynamo has benefited from 
important progresses during the last decade, thanks to 
the combined breakthroughs of observational and modeling 
techniques. Helioseismological studies are now able to re- 
veal the internal velocity field of the Sun (Schou et al. 1998) 
and recent MHD simulations give credence to the idea that 
the dynamo processes at the origin of the large-scale so- 
lar field are mainly operating within a thin interface layer 
separating the radiative core from the convective envelope 
of the Sun : the tachocline. This region allows an efficient 
transformation of a seed poloidal field into a toroidal com- 
ponent as well as a storage of the field on a timescale of the 



order of the solar magnetic cycle (Rempel et al. 2000) . The 
poloidal component of the field is believed to be regenerated 
(with opposite polarity) through the action of the Coriolis 
force, by means of a still mostly unknown mechanism called 
"a-effect", which may also be mostly efficient within the 
tachocline (Dikpati & Oilman 2001). 

The magnetic field observed at the photospheric level of 
the Sun was first detected in sunspots (Hale 1908) which are 
the largest magnetic structures of the solar surface (though 
only covering about 10~* of the photosphere at solar max- 
imum), featuring mostly radially oriented field lines. Mag- 
netic elements are now observed at much smaller scales in 
the so-called quiet photosphere (see Solanki 2001 for a re- 
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view of the properties of such magnetic regions). The in- 
terplay, during the solar cycle, between the small magnetic 
structures and the large-scale field has not yet been under- 
stood, whereas small magnetic elements store a significant 
part of the magnetic energy of the photosphere, and despite 
the fact the distribution of some of these structures has a 
temporal evolution correlated to the solar cycle. Some stud- 
ies also suggest that a turbulent dynamo, disconnected from 
the large-scale solar dynamo, may efficiently generate small- 
scale magnetic structures (Cattaneo 1999). 

Surprisingly, these properties of the solar photospheric 
field are far from the general features observed on other 
active stars, at least among the fast rotator group. The sur- 
face magnetic field of several fast rotating stars has been 
mapped by means of Zeeman-Doppler Imaging (thereafter 
ZDI, Semel 1989). For these extremely active objects, a field 
is generally observed at basically all locations of the photo- 
sphere (see e.g. Donati et al. 2003a), in the form of discrete 
structures that present an homogeneous orientation of field 
lines and cover a significant fraction of the stellar surface. On 
some objects, the surface field is dominated by an azimuthal 
component distributed in arcs or even in complete rings en- 
circling the rotation ajcis at different latitudes. These ob- 
servations suggest that the dynamo operating within these 
objects is very different from that of the Sun, since the az- 
imuthal component of the solar field is supposed to be deeply 
buried at the base of the convective envelope. Its presence 
at the photospheric level on fast rotators therefore suggests 
that their dynamo may be active very close to the surface. 

The observing effort engaged during the last decade to 
investigate, by means of high-resolution spectropolarimetry, 
the magnetic activity of active late-type fast rotators has 
been concentrated on several pre-main sequence stars and 
on the evolved close binary system HR 1099. The present 
study focuses on HD 199178, a member of the small FK Com 
group. Objects belonging to this class are single fast rotating 
late-type giants displaying signs of a strong magnetic activ- 
ity. Their short rotation periods (of order of a few days) 
suggest that FK Com stars have undergone very unusual 
mechanisms during their evolution. In particular, it was pro- 
posed by Bopp & Stencel (1981) that they may result from 
the coalescence of contact binaries. Long-term photometric 
studies of FK Com itself (Jetsu et al. 1991, Korhonen et al. 
2002) report that its cool spot distribution is asymmetric 
(a so-called "active longitude" dominating the spot activ- 
ity), with occasional 180° shifts of the active longitude (a 
behavior dubbed "flip-flop") resulting in a 6.5 years cycle. 
The distribution of cool spots of HD 199178 has also been 
investigated by means of Doppler imaging (Strassmeier et 
al. 1999, Hackman et al. 2001). 

The aim of the present study is to enrich our knowledge 
of this object by reconstructing its magnetic topology alto- 
gether with its spot distribution, in a twofold purpose. We 
first plan to compare the general properties of the photo- 
spheric magnetic field of this evolved star with the magnetic 
topologies mapped for younger objects. We also intend to 
check whether active longitudes reported by Hackman et al. 
(2001) are also observed in our own maps of the cool spot 
distribution and test whether they have a counterpart in the 
magnetic images. 

We first summarize the series of spectropolarimetric ob- 
servations used for this study and describe the modeling 



procedures employed to reconstruct the magnetic topology 
of the star at every observing epoch. We then present the 
main characteristics of its surface brightness and magnetic 
structures and study their evolution on different timescales 
(from a few weeks to several years) , including an analysis of 
the surface differential rotation. We finally summarize the 
main results and outline the informations they provide on 
the nature of the underlying dynamo processes. 



2 OBSERVATIONS, DATA REDUCTION AND 
IMAGING PROCEDURE 

2.1 Observations 

All the spectropolarimetric data presented in this article 
were obtained with the MuSiCoS spectrograph (Baudrand 
& Bohm 1992) fiber fed by its Cassegrain-mounted polari- 
metric module (Donati et al. 1999) . The data reduction, per- 
formed with ESpRIT (Donati et al. 1997), is similar to that 
described by Petit et al. (2003) for HR 1099, including the 
additional wavelength calibration using telluric lines. The 
data sets (listed in Tab. and |5J were obtained between 
1998 Dec. and 2003 Aug., yielding a total number of 380 
brightness spectra and 94 circularly polarized spectra. 

Least-Square Deconvolution (thereafter LSD, Donati et 
al. 1997) was employed to perform a simultaneous extraction 
of the signal from all photospheric spectral features of the 
echelograms. A line mask corresponding to a G5 spectral 
type (Strassmeier et al. 1999) yielded a multiplex gain of the 
order of 30 for the Stokes V profiles. The S/N of LSD Stokes 
I profiles is limited to about 1100 at best, indicating that 
the convolution model underlying LSD is not adapted above 
this accuracy level. Depending on data quality, the multiplex 
gain undergoes small fluctuations, so that raw spectra of 
equal S/N can produce LSD profiles with slightly different 
noise levels. 

The data quality of different observing periods is very 
uneven (due to weather conditions), so that the S/N of LSD 
profiles ranges from about 3,600 (in average in 2001 Jul.) 
to more than 7,000 on 2003 Jul. 31 (in optimal seeing and 
transparency conditions). The large data sets of summers 
2002 and 2003 are split into several subsets (separated by 
blank lines in Tab. ^andj^J. Each subset is used to re- 
construct an individual brightness and magnetic image of 
the star (Sect. 12.21 1. ensuring that each map is computed 
from observations obtained over a time-span consistent with 
the typical short-term evolution of photospheric structures 
(Sect.ESJ. 

2.2 Imaging procedure 

All magnetic images described hereafter are obtained with 
the ZDI code developed by Brown et al. (1991) and Do- 
nati & Brown (1997), following the maximum entropy image 
reconstruction algorithm of Skilling & Bryan (1984). This 
imaging procedure was tested for various stellar parameters 
and observing conditions by Donati & Brown (1997), from 
a series of numerical simulations. They demonstrated that 
the orientation of field lines within magnetic regions can 
be accurately reconstructed for noise levels similar to that 
available in the present study and for inclination angles of 
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Table 1. Journal of observations from 1998 Dec. to 2002 Jul.. Each line corresponds to a full polarization cycle. Columns 2 and 3 list the 
date and hour of observation (first and last exposure of the cycle). Column 4 contains the number of unpolarized/polarized exposures. 
Column 5 lists the total exposure time of each Stokes I individual sub-exposure. We also list the S/N ratios (per 4 kms"'^ velocity bins) 
of the unpolarized and polarized spectra (in columns 6 and 8 respectively, with minimum and maximum values in the sequence) and 
in the associated mean LSD profiles (columns 7 and 9). The multiplex gain between the raw polarized spectra and the mean Stokes V 
profiles is reported in the last column. 
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Table 2. Same as Table [Tl for summer 2003 observations. 
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the rotation axis similar to that derived later in this section 
for HD 199178. However, in case of images reconstructed 
with data presenting an incomplete phase sampling of the 
star, only a partial reconstruction of the magnetic field is 
achieved, containing radial/meridional field regions closest 
to the observed longitudes and azimuthal field structures 
located about 0.1 rotation cycle away from the observed 
longitudes. With sparse phase coverage, the reconstructed 
position of active regions located next to phase gaps can 
also be slightly shifted toward higher latitudes and unob- 
served phases (Petit et al. 2002). Comparisons between ZDI 
images of HR 1099 reported by Donati et al. (2003a) and 
Petit et al. (2003) demonstrate however that maps of the 
same object obtained with simultaneous observations, but 



with different instrumental setups, are very consistent de- 
spite very different noise levels and phase sampling. 

We model the photospheric brightness inhomogeneities 
by means of the two component description of Cameron et 
al. (1992). The average intrinsic profile used for computing 
brightness images is a synthetic Gaussian line reproducing 
the characteristics of a MuSiCoS LSD Stokes I profile of the 
KG star (3 Gem. This option was adopted according to the 
results of Unruh & Cameron (1995), who demonstrated that 
Doppler images reconstructed from a Gaussian line were al- 
most indistinguishable from that obtained using a standard 
star. The template profile was scaled by a factor 0.5 and 1, 
for the spotted areas and the quiet photosphere respectively, 
to mimic the different LSD signatures expected from these 
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two regions (Donati & Cameron 1997). The adopted tem- 
peratures are 5500 K and 4000 K, respectively for the quiet 
photosphere and the spotted regions, both values being close 
to that determined by O'Neal et al. (1998) from TiO band 
modeling. 

The large volume of data we collected for the present 
study allows us to obtain estimates of the imaging parame- 
ters associated with HD 199178. To estimate the radial ve- 
locity Drad, we recoustruct several images from the same 
data set, tuning the value of «rad and choosing the one that 
minimizes the information content of the image (following 
Cameron & Unruh 1994). In order to estimate the projected 
rotational velocity usini, we prefer to adopt the more reli- 
able estimate provided by a classical fit to the data (Donati 
et al. 2003a), from which we can obtain a better adjustment 
of the profiles wings, at the cost of a slightly higher spot 
(or magnetic field) coverage. The values we obtain for -Urad 
and vsini are respectively equal to —28.5 ± 0.5 kms~^ and 
70 ± 1 kms~^, in good agreement with previous estimates 
of Strassmeier et al. (1999). 

The rotation period (or the rotation law in the case of a 
differentially rotating surface) can be estimated in a similar 
way as Vrad, but we leave the detailed analysis of surface 
differential rotation for Sect. 3] AH the images presented 
hereafter are computed assuming that the surface under- 
goes differential rotation during data collection. We assume 
in the reconstruction process a surface rotation law of the 
type: 

n(0 = ficq - dHsin^/ (1) 

where is the rotation rate at latitude /, f2eq the rotation 
rate of the equator (set to 1.934 radd~^, see Sect. and 
dn the difference in rotation rate between the pole and the 
equator (set to 66 mradd"^ by default). We assume a mean 
rotation period equal to 3.3 d (1.904 radd~^) for calculating 
the rotational phases. Given the differential rotation param- 
eters assumed for the star, this period is the one we expect 
at latitude 42° (rotation periods of the stellar surface vary- 
ing from 3.25 d at the equator to 3.36 d in the polar region). 
The adopted ephemeris is : 

JD = 2451150.675 + 3.3 <^ (2) 

where (p is the rotational phase and JD the Julian date. 
Considering this 3.3 ± 0.05 d rotation period and using the 
Dsini value derived above, we estimate the stellar radius to 
be i?.sini — (4.5 ± 0.1) Rq. This value is compatible with 
a previous estimate of Hackman et al. (2001, based on Hip- 
parcos parallax measurements) that yields a stellar radius 
of 4 _R0 < i? < 5 Rq). However, high values of i are obvi- 
ously required to stay consistent with the Hipparcos radius. 
Considering several Doppler images reconstructed with dif- 
ferent values of the inclination, the highest possible value of 
i is about 50°, since higher inclination angles produce ob- 
viously spurious structures on the images. For this inclina- 
tion the value of the stellar radius is roughly 5.8 Rq, higher 
than that proposed by Hackman et al.. However, as empha- 
sized by these authors and by Strassmeier et al. (1999), the 
Hipparcos-based estimate critically depends on the physi- 
cal models used to describe such a post-main sequence star. 
Given the rather unusual physical properties of HD 199178 
(among which its high rotation rate and its uncertain evolu- 
tion on the main sequence, since this single star may come 



from the coalescence of a close binary system), this kind of 
discrepancy is not surprising, therefore we choose to use our 
estimate i = 50°. 

To summarize, stellar parameters adopted in the imag- 
ing process are i — 50°, t^rad = —28.5 kms~^, t;sini= 
70 kms"\ J7eq= 1.934 radd"^ and dn= 66 mradd~\ 



3 RECONSTRUCTED IMAGES 

Nine brightness and eight magnetic images were recon- 
structed from the data sets described in Sect.|5| correspond- 
ing to epochs 1999.02, 2001.54, 2001.97, 2002.46, 2002.50, 
2002.56, 2003.51, 2003.55 and 2003.58 (Figures[Tlto|n|l. The 
early observations of 1998 December provided the first de- 
tection of a magnetic field on HD 199178. Several images 
were reconstructed from subsets of the large data sets se- 
cured in summers 2002 and 2003, to take into account the 
short-term variability of the photosphere occurring during 
data collection (see a discussion of this point in Sect 13.31 . 

Owing to the relatively low S/N of the 2001.54 data set 
(3600 in average in Stokes V, i.e. 40% below the mean S/N 
level of epoch 1999.02 for instance), none of the polarized 
profiles secured in this epoch provide a magnetic field de- 
tection, thus only a brightness image was constructed (Fig. 
I^J. All the surface topologies described hereafter are recon- 
structed with a reduced (hereafter x? ) of 0.65 and 0.9 for 
brightness and magnetic images respectively. A Xr smaller 
than unity is adopted to take into account the fact that 
error bars derived from MuSiCoS observations are slightly 
overestimated (Wade et al. 2000). 



3.1 Brightness surface structures 

The most obvious characteristic of the brightness topology 
of HD 199178 is the presence of a large polar spot in all 
images. This cool region is always quasi-axisymmetric and 
covers most of the stellar photosphere above latitude 60°. 
At several epochs, the polar spot is partly fragmented, with 
smaller spots clearly distinguishable from the main polar 
component (e.g. at epoch 1999.02, when a smaller, secondary 
spot was visible at phase 0.0 and latitude 80°). 

In addition to the polar spot, smaller spots appear at 
lower latitude, but only about 5% of the overall spot cover- 
age is contained within the low latitude features (below lat- 
itude 30°). The latitudinal spot distribution is rather stable 
from one epoch to the next (see Fig. 1101 for epoch 2002.50) 
and repeatedly presents a peak of spot occupancy between 
latitudes 15° and 30°, with a relative gap of spottedness in 
regions next to the equator and within a band extending 
from latitude 30°to the outer limit of the large polar cap. 

In Fig. 1111 we plot the latitudinally-averaged spot cov- 
erage, as a function of the rotational phase, for each subset 
of summer 2002 (left panel) and summer 2003 (right panel). 
The first information provided by these curves is that the 
observed spot distribution is highly varying between near- 
in-time subsets. Part of this apparent evolution can be ex- 
plained by the different phase sampling achieved for different 
subsets, producing a lack of signal at some longitudes. How- 
ever, if phase gaps can be a problem at some epochs (e.g. 
2002.46), most of the subsets possess good phase sampling. 
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Figure 1. Reconstructed images of HD 199178 at cpocli 1999.02, in flattened polar view. The concentric circles correspond (starting 
from the outside) to parallels of latitude —30°, 0° (equator, bold line), +30°ajid +60°. The upper- left panel corresponds to a brightness 
image, while the three other charts show the components of the magnetic field (in Gauss) in spherical coordinates, i.e. radial, azimuthal 
and meridional components of the field in the upper-right, lower-left and lower-right charts respectively. Stokes I and V normalized 
profiles are plotted in the lower part of the figure (left-hand panel and right-hand panel, respectively). Thin lines represent the observed 
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Figure 3. Same as FigEJfor the 2001.97 data set. 
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Figure 5. Same as Fig for the 2002.50 data set. 
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Figure 6. Same as Fig for the 2002.56 data set. 
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Figure 7. Same as FigEJfor the 2003.51 data set. 
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Figure 8. Same as Fig for the 2003.55 data set. 
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Figure 9. Same as Fig for the 2003.58 data set. 
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Spot distribution, 2002 Jul. 



Spot distribution, 2003 Jul. 




Rotational phase Rotational phase 

Figure 11. Distribution of the fractional spot coverage as a function of the rotational phase in 2002 and 2003 (left-hand and right-hand 
panels, respectively). Individual curves are derived from successive data subsets (2002.46 and 2003.51 in full line, 2002.50 and 2003.55 
in dashed line, 2002.56 and 2003.58 in dot-dashes). 
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Figure 12. Same as Fig llll for the fractional magnetic flux. 



particularly in 2003. In order to obtain a more precise esti- 
mate of potential biases due to phase gaps, we compare suc- 
cessive subsets secured at close-by epochs with a different 
phase sampling (but bearing in mind however that short- 
term evolution of the spot distribution also accounts for a 
part of the observed differences, see Sect. 13.31 . 

The spot coverage was marginally higher in 2002 on 
phases ranging from 0.25 to 0.65 (with two distinct maxima 
of spot coverage visible around phases 0.35 and 0.55). In 
2003, the maximum spottedness is rather concentrated be- 
tween phases 0.8 and 0.2. In 2002, longitudes with the largest 
filling factor are clearly linked to the presence of low-latitude 



spots (see Fig.3]to|HJ, while in 2003 features spatially asso- 
ciated to the large polar cap are mostly responsible for the 
higher spottedness around phase 0.0. For other epochs, the 
biases due to uneven phase coverage cannot be as accurately 
estimated as for the 2002 and 2003 data sets. However, if the 
phase sampling is far from optimal at epochs 1999.02 and 
2001.54, there is no significant phase gap at epoch 2001.97, 
for which we note a maximum spot coverage around phase 
0.9 (corresponding plot not shown here). 
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Figure 2. Same as FigQ]for the 2001.54 data set. There was no 
Zeeman signature detected at this epoch, owing to a high noise 
level of Stokes V profiles (see Table 111 . 

3.2 Magnetic topology 

The magnetic topology of HD 199178 features several re- 
gions of intense field strength (in excess of 1 kG). Owing 
to the fact that polarized profiles used here suffer from a 
relatively high noise level (compared to previous studies of 
active stars based on ZDI), the constraints provided on the 
magnetic maps are not very high, so that several assump- 
tions on the field are compatible with an image reconstruc- 
tion at a fixed x? (equal to 0.9). For instance, the maps 
can be calculated assuming either a purely radial field or 
a mixture of all three possible components of the field (ra- 
dial, azimuthal and meridional). To select between these two 
possibilities, the very principle of maximum entropy imag- 



Cool spot distribution 




20 40 60 80 
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Figure 10. Distribution of the fractional spot coverage as a func- 
tion of latitude. The full line curves represent the flux averaged 
over all observing periods. The two dotted curves surrounding the 
mean curve show the standard deviation. 



ing tells us that the image exhibiting the lowest information 
content is the most likely (the information contained in the 
map being calculated as the integrated field strength over 
the stellar surface, thereafter Bmod). Taking epoch 2003.58 
as an example, a purely radial photospheric field leads to 
Binod=74 G. For the same data set, Bmod is equal to 46 G for 
an image containing all three components of the field. Very 
similar results can be derived for all other epochs, with Bmod 
about 30 to 60% higher for a purely radial field than for the 
map containing all three field components. At this stage, it 
can be argued that assuming the presence of an azimuthal 
component is not absolutely necessary for computing the 
magnetic maps. However, several additional reasons lead us 
to consider this component as genuine, as detailed below and 
in Sect. 

The reconstructed azimuthal field component can it- 
self be divided into a couple of subcomponents of opposite 
polarities, with a predominance of clockwise field at high 
latitude, the lower latitudes being dominated by a counter- 
clockwise magnetic field. The left panel of Fig. ll4l shows the 
longitudinally-integrated distribution of the azimuthal field 
(averaged over all observing epochs). Both sub-components 
can be clearly distinguished, with the field polarity chang- 
ing sign around latitude 50°. The fact that the azimuthal 
component presents a consistent large-scale geometry be- 
tween different images (despite very different S/N or phase 
coverage) suggests that the observed structures cannot be 
attributed to a problem of data quality (like noise or phase 
covering) . 

The clockwise sub-component shows up in all images. 
At epoch 2001.97 for instance, a succession of several blobs 
of intense field strength draws the outline of a ring encir- 
cling the pole at the outer border of the cool polar spot. 
We also see an almost complete ring at epochs 2003.55 and 
2003.58, whereas images reconstructed at other observing 
periods feature only a few spots of clockwise field, most of 
the time coincident with the edge of the polar spot. At some 
epochs (e.g. 2003.58), it can be noticed that spots belonging 
to the clockwise sub-component show up within the dark po- 
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lar cap, whereas we expect very few polarization signal from 
this region because of the important brightness contrast with 
the unspotted photosphere (at a fixed field strength, Zeeman 
signatures formed inside the polar cap will be much fainter 
than that formed in the surrounding photosphere). Some 
of these structures may indeed belong to the outer region 
of the polar spot, but reconstruction biases of the imaging 
code can slightly shift these magnetic regions toward higher 
latitudes (Petit et al. 2002). It cannot be excluded how- 
ever that Zeeman signatures are observed inside the polar 
spot. In this case, the reconstructed magnetic field must be 
considered with a caution, since the field strength is likely 
underestimated. 

The counter-clockwise component of the azimuthal 
magnetic field appears more discreetly between latitudes 10° 
and 40°. At epoch 2003.58 for instance, several spots of az- 
imuthal field appear around latitude 40°, but with a rela- 
tively low field strength (reaching about 0.6 kG at phases 
0.5 and 0.9). The low level of the magnetic flux emerging 
from these regions (Fig. I14II explains why such components 
cannot (or only marginally) be reconstructed at epochs for 
which the data suffer from a lower S/N (epochs 2002.50 and 
2002.56 for instance). 

In addition to areas hosting an azimuthal magnetic 
field, we reconstruct several regions in which the field is 
mostly radially oriented. On most images, the radial field 
is organized in an intricate pattern mixing opposite polari- 
ties, with a large-scale structure much less obvious than that 
of the azimuthal component. We note also that the radial 
field regions are varying a lot between images corresponding 
to close-by epochs, which may lead to question their reality. 
Numerical simulations suggest indeed that the noise pattern 
is preferentially reconstructed as radial field regions (but this 
effect is expected to be significant only for S/N much lower 
than that achieved here). As a test, we divide all the epoch 
2003.58 data into two independent data sets by taking the 
odd numbered spectra to generate one image and the even 
numbered spectra to generate a second image (both images 
not shown here). The smallest radial field regions obtained 
in the corresponding images show little similarities between 
both maps, while all larger structures are coherently recon- 
structed. We emphasize however that reconstructing a map 
with one half of the profiles is tantamount to increasing the 
noise by a factor \/2. Given the fact that the reconstructed 
magnetic distribution is therefore less accurate, we caution 
that discrepancies at small scales in the radial field pattern 
do not demonstrate that small-scale structures are not re- 
liable in images better constrained. Furthermore, the good 
consistency achieved for all large structures is a good indi- 
cator of their reality. 

The largest radial field regions are observed at epochs 
2002.46 and 2003.58, when two large unipolar areas are seen 
close to the pole. Contrary to the azimuthal field, no lati- 
tudinal organization of the radial field can be readily seen 
on the different images. However, if we average the radial 
field component over all observing epochs (Fig. 1141 , we note 
that a positive field dominates at high latitude, while re- 
gions presenting a negative polarity are mostly confined be- 
low latitude 40°. Despite this statistical trend, it can be 
seen on some images that a negative polarity can sometimes 
dominate the radial component at high latitude (at epochs 
1999.02 and 2002.56 for instance). 



A small part of the polar magnetic field is sometimes 
reconstructed as a meridional field component (see e.g. Fig 
ISJ. These high latitude azimuthal and meridional features 
indeed belong to the same physical structure (a region con- 
taining mostly horizontal field and located over the pole), an 
effect (detailed by Donati 1999) that can be attributed to the 
singular nature of the pole in spherical coordinates. Other 
structures appearing in the meridional field component at 
lower latitude (see, e.g., very faint structures around phase 
0.5 at epoch 2003.58) can be partly produced by known re- 
construction cross-talks between radial and meridional fields 
(Donati & Brown 1997), an effect essentially negligible in the 
case of HD 199178 thanks to its relatively high inclination 
angle. 

The radial magnetic field is always the main contributor 
to the photospheric magnetic energy (at roughly 85% of the 
overall energy). However, numerical simulations suggest that 
this fraction is sensitive to the S/N and may become high 
for low-quality data (part of the noise producing radial field 
structures, even for very slight over-fitting of the data). In 
order to avoid this problem, a solution consists in calculating 
the energy stored in the longitudinally-integrated field. By 
doing so, we average out most of the noise, but we also 
ignore most of the energy contained in small-scale magnetic 
regions. It is however an opportunity to estimate the energy 
stored within the large-scale axisymmetric field, which is 
otherwise partly hidden among small magnetic features. We 
find that the radial field contributes only to some 10% of the 
global energy contained in the large-scale field (at epochs 
2001.97 and following), which means in other words that 
the axi-symmetric field is largely dominated by its azimuthal 
component. 

At epochs 2002.46, 2002.50 and 2002.56, the low- 
latitude ring of azimuthal field steadily contains about 14% 
of the total amount of energy stored within the azimuthal 
component of the axisymmetric field. This fraction is stable 
despite large variations of the S/N and phase coverage over 
the 3 corresponding data sets, suggesting that data qual- 
ity only marginally impacts this estimate. The fraction was 
equal to 6% six months earlier at epoch 2001.97 (with signif- 
icantly better S/N and comparable level of phase coverage) 
and reaches about 20% between 2003.51 and 2003.58. We 
suggest that this evolution may be genuine, with respective 
weights of both azimuthal field rings remaining rather sta- 
ble over a few weeks, but showing significant evolution on 
timescales ranging from a few months to several years. 

The latitudinally-averaged distribution of the photo- 
spheric field as a function of the rotational phase is plotted 
in Fig. 1121 for data sets secured during summers 2002 and 
2003. A larger magnetic fiux is observed between phases 0.2 
and 0.4 during summer 2002, but with significant differences 
between the three curves corresponding to epochs 2002.46, 
2002.50 and 2002.56. We wiU discuss in Sect. lO whether 
short-term evolution of active regions can account for this 
apparent discrepancy. At epoch 2002.56 and during summer 
2003, the distribution of the magnetic flux is mostly inde- 
pendent of the rotational phase. In particular there is no 
apparent correlation between the phase distributions of cool 
spots (Fig. im and magnetic field. 
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Figure 13. Stokes I profiles of HD 199178 obtained at phase 
0.5588 on 2002 June 13 (dashed fine), and at phase 0.5629 on 
July 26 (full Une). 



3.3 Rapid surface changes 

During summers 2002 and 2003, HD 199178 was observed 
every clear night, thus providing a monitoring (over two dis- 
tinct periods covering about 45 d) of the evolution of surface 
structures. 

The intrinsic evolution of the photosphere is first obvi- 
ous in the profiles themselves. For instance, striking differ- 
ences are visible between the brightness profiles secured on 
2002 June 13 and 2002 July 26, i.e. 43 days apart (Fig.lTa. 
despite a difference in rotational phase as small as 0.4% of a 
rotation cycle. In this specific case, the wide bump appear- 
ing at the center of the most recent profile has no equivalent 
6 weeks before. 

The second evidence for surface variability arises from 
a monitoring of the Xr of the reconstructed images. When 
grouping the data secured from 2002 June 28 to 2002 July 
28 in a single data set (i.e. grouping epochs 2002.50 and 
2002.56), we cannot reconstruct surface images with a Xr 
lower than 0.85 and 1 for the brightness and magnetic im- 
ages respectively, as opposed to a Xi oi 0.65 and 0.9 when 
reconstructing separate images from each sub-set. To ensure 
that the data sets employed to produce the surface images 
do not spread on too long a timespan, we split the global 
data sets of summers 2002 and 2003 into successive subsets 
corresponding to epochs 2002.46, 2002.50, 2002.56, 2003.51, 
2003.55 and 2003.58 (respectively covering 14, 13, 11, 15, 13 
and 11 nights). We then use this time series to get a direct 
view at the surface evolution of brightness and magnetic 
structures. In the following paragraphs, we voluntarily limit 
the comparisons between pairs of images to high-latitude 
structures and to the low-latitude regions located at rota- 
tional phases observed in both data sets, in an attempt to 
avoid as much as possible difi'erences produced by recon- 
struction biases associated to phase gaps. 

Considering the set of brightness images, we first note 
the continuous birth and disappearance of low latitude 
spots. For instance, the spot located at phase 0.4 and lat- 
itude 70° at epoch 2003.58 cannot be associated to any 
structure observed at epochs 2003.51 and 2003.55. On the 
other hand, the spot reconstructed around phase 0.17 and 
latitude 25° at epochs 2003.51 remains clearly visible until 



2003.58 (small changes in its exact location and shape from 
one epoch to the next being compatible with usual recon- 
struction biases). 

Rapid changes also afi^ect brightness inhomogeneities at 
the limit of the polar spot during similar timescales. The 
edges of the polar cap (from phases 0.0 to 0.5) are rapidly 
evolving during the summer of 2002, with two secondary 
spots appearing next to the main component of the polar 
spot between epochs 2002.46 and 2002.50, then followed by 
a third region, located around phase 0.98 at epoch 2002.56. 
Note that these small spots are located around latitude 60° 
and therefore are sufficiently far from the pole to produce 
significant rotational modulation in the spectra. Considering 
also that these structures are never eclipsed during stellar 
rotation, we can conclude that the location of such cool spots 
is very well constrained, so that the evolution we report here 
is very likely to be real. 

The same type of evolution is observed for magnetic re- 
gions, with striking changes in the magnetic topology occur- 
ring on timescales as short as a couple of weeks. One of the 
most obvious change concerns the large azimuthal field re- 
gion located at phase 0.5 and latitude 45° at epochs 2003.55 
and 2003.58, which has no counterpart at epoch 2003.51. 
Another striking example concerns the very intense, large 
spot located at phase 0.4 and latitude 60° at epoch 2002.46. 
The field strength inside this region, in excess of 1 kC at 
epoch 2002.46, rapidly decreases to about 100 G in epoch 
2002.56 (this evolution explains part of the differences be- 
tween curves plotted in the left panel of Fig. I12II . We cau- 
tion however that the relatively low S/N of data sets corre- 
sponding to epochs 2002.50 and 2002.56 may also be partly 
responsible for the observed evolution. 

Concerning the spots of radially-oriented magnetic field, 
only the largest structures seem to possess lifetimes larger 
than a couple of weeks. As an illustration, the large region of 
positive field located near the pole at epoch 2002.46, which 
has already partly disappeared at epoch 2002.50, is then 
replaced by a region of opposite polarity at epoch 2002.56. 
We must nevertheless keep in mind that the location and 
field strength of magnetic regions inside the dark polar spot 
are likely to be affected by reconstructions biases, as already 
outlined in Sect. 13.21 



4 DIFFERENTIAL ROTATION 

In this section, we study the motion of surface structures un- 
der the effect of differential rotation, applying the dedicated 
method of Petit et al. (2002) and assuming a surface rotation 
law similar to that given by equation Considering a time- 
series of profiles, up to 2,000 images are computed (either 
brightness or magnetic maps, depending whether polarized 
or unpolarized profiles are considered) assuming different 
values of the differential rotation parameters f2cq and dQ 
and imposing a constant information content in all images. 
In the case where the data quality is good enough, the re- 
constructed images present a x? minimum corresponding to 
the most likely set of differential rotation parameters. Con- 
sidering a Xi map obtained from the whole set of computed 
images (e.g. Fig. I15L the region around the x? minimum 
can be approximated by a paraboloid, which curvature ra- 
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Figure 14. Latitudinal distribution of tlie azimuthal (loft) and radial (right) components of HD 199178 magnetic flux. The full line 
curves represent the flux averaged over all observing periods. The dotted curves show the standard deviation. 



Table 3. Surface differential rotation parameters derived for 
HD 199178, using Stokes V profiles. For each epoch, f2oq (equa- 
torial rotation rate) and dQ (difference of rotation rate between 
equator and pole) are listed. 

Date f2eq dQ 

year radd~^ mradd"^ 

2001.97 1.934 ±0.018 66 ± 28 
2002.48 1.932 ± 0.033 79 ± 61 
2003.57 1.931 ±0.017 86 ± 20 



dius around its minimum gives the formal error bars on both 
parameters Qcq and dQ (Donati et ah 2003b). 

Only three of our data sets provided a detection of dif- 
ferential rotation (Table I^J. The first one corresponds to 
observations obtained at epoch 2001.97. The two others are 
subsets taken from the large data sets of summers 2002 and 
2003. In summer 2002, we consider profiles secured between 
June 15 and July 2. In summer 2003, the subset is built 
from observations ranging from July 21 to August 5. Sub- 
sets of 2002 and 2003 observations are chosen to present the 
best compromise between phase sampling, noise level and 
time-span. In particular, we ensure that the time-length of 
all data sets used to derive the differential rotation parame- 
ters is consistent with the typical timescale of local changes 
in the photospheric structures, as discussed in Sect. 13.31 
Measurements of the dQ parameter give similar results in 
2001.97 and 2003.57, with dQ equal to 66 ± 28 mradd"^ 
and 86±20 mradd"^ respectively, indicating the detection of 
a solar-like differential rotation (the equator rotating faster 
than the pole) to within 4.3cr. The third estimate, corre- 
sponding to epoch 2002.48, suffers from error bars 2 to 3 
times larger, but stays consistent with measurements from 
the two other epochs. 

No other data set constituted of Stokes V profiles could 
provide a detection of differential rotation (no paraboloid 



shape of the Xr map around x? minimum). A higher noise 
level, a sparser phase sampling or a local evolution of the 
distribution of magnetic regions during data collection may 
be responsible for this failure. We also note that no detec- 
tion was achieved using Stokes I profiles. This is not a sur- 
prise, since cool spots are mostly concentrated close to the 
pole on HD 199178, few tracers being available at lower lat- 
itude. On the contrary, the magnetic topology offers tracers 
densely distributed over the stellar surface, therefore span- 
ning a larger range of rotation periods and allowing an easier 
detection of a surface shear. 

As outlined by Petit et al. (2002 and 2003), measure- 
ments of differential rotation can be partly affected by biases 
arising from low data quality. Also problematic is the po- 
tential aliasing due to intrinsic evolution of the photosphere 
(a newborn active region can be mistaken for an older one 
that has vanished during data collection, leading to inade- 
quate evaluations of the shear). However, it is highly unlikely 
that an aliasing problem occurs for the three data sets pre- 
sented here, and even more unlikely that the produced bias 
could shift our estimate of differential rotation toward the 
same (spurious) value. Given the fact that the data sets also 
present very different S/N and phase sampling, the three 
consistent measurements presented here suggest that the ob- 
served shear is genuine. 



5 DISCUSSION 

5.1 Brightness distribution 

The magnetic activity of HD 199178 has been monitored 
by the present study over 4.5 yr. Several characteristics of 
its brightness topology remain constant during this whole 
period, the most obvious being a large axisymmetric po- 
lar spot reconstructed above latitude 60° . While dark polar 
caps are known to be a usual feature on Doppler images 



Magnetic field & differential rotation of HD 199178 19 



map, 2003 Jul. 

T I I I I I I I I I I r 




1,9 1,95 
(rad/d) 

Figure 15. reduced map in tlie differential rotation param- 
eter plane, obtained from polarized data of epoch 2003.57. Qcq 
is the rotation rate of the equator and dfl the difference in rota- 
tion rate between the pole and the equator. The Icr limit on the 
two parameters (considered separately) corresponds to the black 
region. 

of fast rotators, their exact shape and evolution with time 
can significantly vary for different objects. The only other 
evolved star for which ZDI has been performed up to now 
(the Kl sub-giant of the RS CVn system HR 1099, see e.g. 
Petit et al. 2003) also possesses a cool polar region, but with 
a more complex shape and a fluctuating location of its cen- 
troid. Taking into account the regular, almost axi-symmetric 
shape of the polar spot observed on HD 199178, Hackman 
et al. (2001) argue that the flat bottom of the photospheric 
Stokes I profiles (interpreted by imaging codes as a polar 
spot) may indeed be partly due to a distortion of line profiles 
under the effect of anti-solar differential rotation (the polar 
region rotating faster than the equator). Our own measure- 
ments of the surface shear of HD 199178 suggest a solar-like 
differential rotation of magnetic structures, arguing against 
this possibility. Moreover, numerical simulations show that 
for shear intensities similar to that we measure, line profile 
distortions do not exceed lO"** of the continuum level, al- 
lowing us to safely assume that the differential rotation of 
HD 199178 does not produce significant biases in the imag- 
ing process. 

5.2 Magnetic topology 

The magnetic topology of HD 199178 presents at all ob- 
serving epochs large regions in which field lines are mostly 
azimuthally oriented. Such magnetic regions have already 
been observed in several objects (e.g. Donati et al. 2003a) 
with sufficiently high S /N to erase any doubt on their reality. 
Here, two long-lived rings displaying opposite polarities are 
observed, encircling the pole at two distinct latitudes. The 



fact we observe this large-scale structure on all maps sug- 
gests that it cannot be attributed to an artifact due to the 
noise pattern. This general configuration of the azimuthal 
field is very similar to that of HR 1099, whereas no such ob- 
vious axi-symmetric structure has been reported for younger 
stars (which suggests that such large-scale characteristics of 
the magnetic topology are not a systematic bias produced 
by ZDI) . The radial field component does not present a sim- 
ilarly axi-symmetric structure. 

This predominant azimuthal component may be con- 
nected to the toroidal component of the large-scale dynamo 
field, by analogy with dynamo models developed for the 
Sun. However, this would imply at the same time that the 
dynamo operating in HD 199178 is possibly different from 
that at work in the Sun. In the solar case, the toroidal field 
is believed to be mostly confined at the interface between 
the radiative core and the convective zone. The observa- 
tion of this toroidal component at the photospheric level 
would therefore suggest that the dynamo of HD 199178 is 
active either very close to its surface, or maybe in its whole 
convective envelope. Steep sub-surface velocity gradients are 
observed in the Sun (Corbard & Thompson 2002) and it was 
investigated by Dikpati et al. (2002) whether a near-surface 
dynamo could be generated in this layer. They concluded 
however that such a mechanism could only marginally con- 
tribute to create a large-scale field. Recent simulations of 
Kiiker & Stix (2001) however suggest that sub-surface veloc- 
ity gradients may be stronger for stars with deep convective 
zones and for fast rotators. In this context, the possibility 
that a sub-surface dynamo may be efficient in a star like 
HD 199178 is an interesting option that could partly ac- 
count for the observed strong azimuthal photospheric field. 

5.3 Active longitudes 

Recent studies (based on long-term photometric observa- 
tions) report the presence of active longitudes on FK Com, 
a star very similar to HD 199178 (Korhonen et al. 2002). 
Furthermore, the same authors report occasional 180° shifts 
of the activity peak longitude, a phenomenon usually called 
"fiip-flop" . The kind of monitoring we propose in the present 
study is not adapted to analyze such an effect, due to an 
inadequate time interval between successive images (a reg- 
ular monitoring of HD 199178 over several months would 
be an appropriate observational basis in this aim). In our 
observations there is only a marginal evidence that some ro- 
tational phases are more spotted than others (Fig. IllH . It 
is also worth noticing that we do not detect any correla- 
tion between the phase distributions of spot coverage and 
magnetic flux (Fig. 1121 . This may suggest that the observed 
photospheric magnetic fleld is not directly connected to the 
internal fleld producing the cool spot pattern. In particular, 
it gives further support to the idea that the observed mag- 
netic structures may be formed close to the stellar surface, 
while cool spots may reveal the action of a dynamo seated 
deeper in the stellar interior. 



5.4 Magnetic cycle 

The observed topology of the magnetic field is very stable 
over the whole observing period. While local modifications 
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of the magnetic pattern occur on tirricscales as short as a 
couple of weeks, some general characteristics of the large- 
scale field (like the rings of azimuthal field) only undergo 
marginal changes on the 4.5 yr observing window. Only two 
noticeable long-term changes can be reported here. 

The first one is the progressive polarity reversal of 
the radial component of the field at high latitude, between 
epochs 2002.46 and 2002.56. No simultaneous variability of 
other large magnetic structures was noticed. A similar phe- 
nomenon was earlier reported for AB Dor (Donati et al. 
2003a) but over a longer timescale, of order of 1 yr. Even if 
we cannot exclude at this stage that this local evolution is 
part of the global activity cycle of HD 199178, the overall 
stability of tfie largest magnetic structures rather suggests 
that this kind of event is only local, thus not necessarily 
connected to a global variability of the large-scale field. The 
second noticeable evolution, which shows up on a timescale 
of several years, is the increasing fraction of magnetic energy 
stored at low latitude in the azimuthal field component. 

5.5 Surface differential rotation 

We report the probable detection of a solar-like 
differential rotation on HD 199178, derived by a moni- 
toring of the relative motion of surface magnetic tracers. 
The intensity of the surface shear, estimated from the 
difference in rotation rate between equatorial and polar 
regions, is roughly solar in magnitude, with a lap>-time 
(time for the equator to lap the polar region by one 
complete cycle) equal to 80 days. The repeated detection of 
differential rotation from independent data sets and the size 
of the related error bars allows us to rule out the possibility 
of an anti-solar differential rotation. 

The fact that the shear intensity is similar in magni- 
tude to that measured on several fast rotating young dwarfs 
(Cameron et al. 2000) may indicate that a difference in the 
depth of the convective zone has only a limited impact on the 
shear level. It was also suggested by Cameron et al. (2000), 
considering a small sample of G and K active stars, that the 
lap-time may decrease for increasing stellar masses. With a 
mass slightly higher than that of other stars of this sample 
(a stellar mass of 1.65M0 was proposed by Ackman et al. 
2001), the lap-time of HD 199178 is also among the shortest 
measured on other active G and K stars, therefore giving 
further support to this idea. 

A very weak solar-like surface shear similar to that de- 
tected on HR 1099 (Petit et al. 2003 report for this star a 
lap-time equal to 480 d) is excluded to the 4cr level. Both 
HR 1099 and HD 199178 are evolved stars with deep con- 
vective envelopes, therefore we do not expect a difference in 
their evolutionary stage to be responsible for such discrep- 
ancy. A possible explanation is that the strong tidal forces 
operating in the convective envelope of HR 1099 are respon- 
sible for its very weak differential rotation, though observa- 
tions of other close binaries are obviously needed to confirm 
this result and disentangle the influence of a tidal torque 
from that of other stellar parameters. 

The lifetime of surface structures, reported to undergo 
fast changes o\'er periods as short as a couple of weeks, is 
similar to values suggested by earlier studies of young dwarfs 
(Barnes et al. 1998), but much shorter than on HR 1099 
where typical active regions remain stable on timescales as 



long as 4 to 6 weeks. The smaller size of most magnetic and 
spotted regions of HD 199178 may partly explain this differ- 
ence in lifespan, as well as its stronger differential rotation 
producing a more intense shear of active regions. 



6 CONCLUSIONS AND PROSPECTIVES 

This study reports the detection of a photospheric magnetic 
field on the FK Com giant HD 199178. The spatial distribu- 
tion of the field is reconstructed by means of ZDI. The mag- 
netic topology of HD 199178 presents at all observing epochs 
large regions in which field lines are mostly azimuthally ori- 
ented. This observation, together with similar ones reported 
for other fast rotators (Donati et al. 2003a), suggests that 
the dynamo processes generating the magnetic activity of 
HD 199178 may be active very close to the stellar surface. 

Local short-term evolution of surface brigfitness and 
magnetic structures are observed on timescales as short as 
a couple of weeks. A polarity reversal is reported during 
summer 2002 for high-latitude magnetic regions hosting a 
radial field. A slow increase of the fraction of magnetic 
energy stored at low-latitude in the large-scale azimuthal 
component is also observed. Beside these local evolutions, 
the largest features of the brightness and magnetic topolo- 
gies remain stable over several years. Although the observed 
surface changes remain mostly limited to small structures 
on our 4.5 yr monitoring, future observations may tell us 
whether the magnetic field of HD 199178 undergoes a cycli- 
cal evolution similar to that observed on the Sun. 

It is finally suggested that the surface of HD 199178 is 
sheared by a solar-like differential rotation. The difference 
in rotation rate between equatorial and polar regions is re- 
ported to be about 1.5 times that of the Sun. 

The observations reported here are at the limit of the ca- 
pacity of the MuSiCoS spectropolarimeter, in terms of S /N 
and spectral resolution. Much more accurate observations 
are therefore expected for the same star when the new gen- 
eration of spectropolarimeters becomes available, with ES- 
PaDOnS at the Canada-France-Hawaii Telescope (Donati 
2003) and NARVAL at TBL (Auriere 2003). Better data may 
in particular allow us to probe more discreet photospheric 
effects, like possible secular fiuctuations in the amount of dif- 
ferential rotation, as already reported for other active stars 
(Donati et al. 2003b). 
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